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Abstract 



We study the multi-exciton optical spectrum of vertically cou- 
pled GaN/AlN quantum dots with a realistic three-dimensional direct- 
diagonalization approach for the description of few-particle Coulomb- 
correlated states. We present a detailed analysis of the fundamental properties 
of few-particle/exciton interactions peculiar of nitride materials. The giant 
intrinsic electric fields and the high electron/hole effective masses give rise 
to different effects compared to GaAs-based quantum dots: intrinsic exciton- 
exciton coupling, non-molecular character of coupled dot exciton wavefunc- 
tion, strong dependence of the oscillator strength on the dot height, large 



ground state energy shift for dots separated by different barriers. Some of 
these effects make GaN/AlN quantum dots interesting candidates in quan- 
tum information processing. 
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I. INTRODUCTION 



In the last few years semiconductor quantum dots (QDs) has been an area of intensive 
research in condensed matter physics [1]. The possibihty of tailoring some of their structural 
parameters represents a fundamental tool in order to study basic physics and to achieve 
technological applications [2] , [3] , [4] . The strength of the lateral or vertical couphng between 
two quantum dots represents a further degree of freedom for engineering optical and electrical 
properties and for studying few-particle phenomena [5] , [6] ; such coupled nanostructures are 
also known as "artificial" or "macro" molecules. 

The growth of GaN/AlN quantum dots [7] and the demonstration of their vertical and 
lateral alignment [8] has been recently achieved. While GaAs-based quantum dots have been 
widely studied both theoretically and experimentally [9], [10], [11], GaN quantum dots are 
becoming a subject of increasing interest for their potential technological applications, such 
as quantum information processing [12], single electron read-out devices [13] or spintronics 
[14]. 

Our analysis is focused on few-particle effects in nitride dots. We address some dis- 
tinguished few-particle phenomena typical of nitride QDs and mainly stemming from the 
built-in giant electric field which characterizes such nanostructures. We analyze the behav- 
ior of the intrinsic exciton-exciton dipole coupling as a function of the various structure 
parameters. In addition our calculations show that the ground state excitonic transition of 
two identical GaN dots is characterized by a strong blueshift when their relative distance is 
decreased. The corresponding ground-state excitonic wavefunction preserves its atomic-like 
behavior, in contrast to what happens in GaAs-based couple QDs, where a redshift of the 
ground-state and a splitting in bonding and anti-bonding molecular-like states is observed 
[15]. We stress that our analysis is also relevant for the experimental realization of the 
quantum information processing strategy proposed in [12], in which a large biexcitonic shift 
is necessary for energy selective addressing of the different few-particle excitation with fs/ps 
laser pulses. 
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The paper is organized as follows. In section II we compare the characteristics of GaAs 
and GaN-based heterostructures, in sec. Ill we develop our theoretical approach; in section 
IV we introduce the investigated system; in section V we focus on exciton-exciton interac- 
tions and on the oscillator strength of the transitions. In section VI we study two coupled 
identical GaN dots and we compare GaAs and GaN macro molecules investigating the role 
of the different parameters involved. In section VII the main results are summarized and 
some conclusions are drawn. 

II. GAN VERSUS GAAS-BASED QUANTUM DOTS 

Let us first of all compare the general characteristics of GaN-based quantum dots with 
the more famihar GaAs-based ones. The GaAS compounds have a zincblende structure 
characterized by a face-centered cubic cell. This structure does not present spontaneous 
polarization, and, in the absence of applied electric fields, the dipole moment of correlated 
electron-hole excitations will be negligible and so will be the interaction between excitons 
created in different QDs. For application purposes however, it would be desirable to have 
a controllable source of interaction between different QDs, since the latter can be naturally 
thought as the different units of a device. A possibility for creating such an interaction has 
been envised and studied in [10,11]. The underlying idea is to create interacting dipoles, 
i.e. to separate electron and hole excitonic components by means of an in-plane externally 
applied (static) electric field. 

The great advantage of III-V nitride compounds, as GaN, is that they may present 
intrinsic electric fields, and, as a consequence, they automatically have built-in such an 
interaction. III-V nitride compounds may present in fact a wurzite type of structure based 
on a hexagonal cell, which is compatible with spontaneous bulk polarization, i.e. they 
present a nonvanishing dipole per unit volume. In the heterostructures we consider, this 
polarization is accumulated at the GaN/AlGaN interfaces due to the asymmetry of GaN 
and AIN unit cells. 
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In the case of QDs a strong strain-induced piezoelectric field along the (0001) direction 
must be added to this effect. This is considerably stronger than in GaAs-based nanostruc- 
tures because the piezoelectric constants in nitrides are orders of magnitude greater than in 
other III-V compounds. The sum of the two contributions results in a strong built-in elec- 
tric field of the order of few MV/cm which is oriented along the growth direction, and has 
opposite sign inside and outside the dot. GaN-based QD are then characterized by strong 
intrinsic electric fields which automatically creates intrinsic dipoles out of each and every 
electron-hole correlated pair excitation. 

Let us consider two vertically stacked QDs coupled by the exciton-exciton interaction 
just described. In GaAs the particle distributions corresponding to the biexcitonic ground 
state generated by laser pulses having the same polarization, will approximately correspond 
to two parallel dipoles (one for every QD) aligned along the (in-plane) field direction [11]. 
In this way excitons in two vertically stacked QDs will have a positive Coulomb interaction. 

If the corresponding GaN based structure instead, such a state is formed by approxi- 
mately two dipoles but stacked in the growth direction. Their interaction energy will then 
be negative. This interaction energy is the so called biexcitonic shift, i.e. the difference in 
energy between the transition energy corresponding to the creation of a certain exciton and 
the one corresponding to the same transition but in the presence of another exciton. 

In GaAs based structures the biexcitonic shift can be tuned by engineering dots param- 
eters as height and base and by modifying the value of the external applied field. In the 
case of GaN based the strength of the built in field is instead determined by the structure 
parameters. In both materials and for experimentally reasonable parameters, it is possible 
to achieve (at least theoretically ) biexcitonic shifts of the order of few meV, which are 
consistent with fast, sub-picosecond, exciton dynamics. 

From the previous analysis it is clear that both GaAs and GaN structures may present 
a tunable exciton-exciton coupling which can be used for different applications and in par- 
ticular to perform two-qubit conditional operations [10-12]. One drawback of applying an 
external electric field (as needed in GaAs based structures) is that a too strong electric field 
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might ionize the trapped charges. An externally applied field also means more complex 
overhead circuits. 

Due to the strength of the built-in electric field, the difference in energy between the two 
lowest excitonic transitions is one order of magnitude bigger in GaN-based structures than 
in the GaAs-based ones. Additionally in the GaN-based structures the oscillator strength 
associated to ground state excitons in QDs of different heights are different because the QDs 
height roughly determines the dipole length and this strongly infiuence electron and hole 
wavef unction overlap. 

III. THEORETICAL MODEL 

The theoretical approach employed to study the optical response of our GaN nanostruc- 
ture is a fully three-dimensional exact-diagonalization scheme, as described in [11]. The 
confinement potential of GaN " macromolecule" is modeled as parabolic in the x — y plane 
and as a sequence of triangular-like potential wells along the growth [z] direction (see Fig.l). 

The physical system under investigation is a gas of electron-hole pairs confined in a 
semiconductor quantum dot or in two vertically coupled quantum dots. The Hamiltonian is 
the sum of three terms: 

H'^ + H'^ + H' (1) 

The term H'^ is a sum of single-particle Hamiltonians which describes a gas of nonin- 
teracting carriers, electrons and holes, and includes the QD quasi-OD confining potential. 
The term H'^'^ describes the correlation of the carriers via the two-body Coulomb interac- 
tion. The term H' is a light-matter interaction Hamiltonian, which accounts for the laser 
light absorption in a quantum dot. We consider the multi-exciton optical spectra, i.e., the 
absorption probability corresponding to the generic N ^ N' transition, where N is the 
exciton number, and in particular we evaluate [11]. the excitonic (0 — > 1) and biexcitonic 
(1 — > 2) optical spectra. The biexcitonic (1 — > 2) optical spectrum describes the creation of 
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a second electron-hole pair in the presence of a previously generated exciton. Here, we shall 
consider parallel-spin configurations only. In the next three subsection we will describe the 
three terms of the Hamiltonian. 



A. Single particle description 

We will work in the usual effective-mass [18] and envelope- function approximation [19]. 
Within such approximation scheme, the non-interacting electron and holes wavefunctions 
are described by the following Schrodinger equation: 



i^e/h + V'^'^^lje/h^Ee/h^e/h (2) 



where e/h describes the set of quantum numbers for electrons (e) and holes (h); is 
the envelope function of the quantum state; ml^^^ is the bulk effective mass for electrons 
(e) and holes (h); E^/h are the energy levels; V'^^'^ is the three-dimensional confinement po- 
tential. Since we are interested only in the lowest energy levels, we can approximate the 
confinement potential as the sum of two potential profiles acting in the parallel and perpen- 
dicular direction to the growth plane. Moreover, a 2D parabolic potential in the x-y plane 
has been proven to reproduce experimental data [9]. Therefore in our model we consider a 
2D parabolic potential, whose analytic solutions are known, and we solve numerically the 
Schrodinger equation along the growth direction using a plane-wave expansion technique 
[20] . We expand the unknown envelope function in a plane- wave basis: 

We use periodic boundary conditions with a box of length L; k — n(27r/L) are the 
reciprocal lattice vectors. By substituting the plane-wave expansion (3) in the Schrodinger 
equation (2) we obtain an eigenvalue equation: 

J2kiHkk' - ESkk')bk' = (4) 
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where Hk}^i are the matrix elements of the single-particle Hamiltonian in a plane- wave ba- 
sis. A direct diagonalization can be done with standard commercial packages. The electron 
and hole states expressed in terms of creation/destruction operators are: 

|e) = 410) 

\h) = 4|0) (5) 

where |0) is the vacuum state. The single-particle Hamiltonian of the system can be 
rewritten as 

H^^H^ + H^ = Y^Eeclce + E^EndU, (6) 



B. Coulomb correlations 

The single-particle eigenfunctions and the corresponding eigenvectors are the input for 
the calculation of Coulomb interactions between the carriers. In writing H"^ we neglect 
Auger processes, that take place only at high particle densities and energy far from the 
bottom of the band. Only processes conserving the number of carriers are considered: 

JJCC JJ^^ -|- JJ^^ _j_ JJ^^ 

2 5Ieie2e3e4^^l«2e3e4cii 42^63^64 + 
2E^lft,fe3^,^'*l'^2fe3fe44l42 4344 + 

~^eie2/ii/i2^^^2'*^'*2'"^i'^'*i'^''2*"^2 ^'^^ 

The physical interpretation of the terms in H"^ is respectively, electron-electron and hole- 
hole repulsive interaction (if^"^ and H^^ ) and electron-hole attractive interaction {H^^), 
while V indicates the the Coulomb potential matrix for a generic two-particle transition. 
After having obtained by direct diagonahzation the 3D single-particle electron and hole 
wavefunctions, we calculate the matrix elements of the complete many-body Hamiltonian 
[H'^ -\- H'^'^ ) in the basis of products of electron and hole eigenstates of the single-particle 
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Hamiltonian. We consider many-body state with the same number of electron and holes, 
that are in general denoted as exciton (N=l), biexciton (N=2), etc. By direct diagonaliza- 
tion we obtain the energies and wavef unctions, that will be expressed as a linear combination 
of products of single-particle states. 



C. Matter- light interaction 

The light-matter interaction Hamiltonian is: 

H' = -E{t) (E,,/^:.44 + Ee./^e.Ce4) (8) 

E(t) is the electromagnetic field of the laser and jieh is the dipole matrix element for the 
e-h transition that can be factorized as a "bulk" and a wavefunction dependent part: 

A^e/i = tJihuik j 'il^l{r)iph{r)dr (9) 

The matrix elements of the light-matter Hamiltonian are different from zero only for 
the transitions N^[N + 1) and {N + 1)^A^ that correspond to the absorption and the 
emission of a photon, respectively (N is the number of electron-hole pairs, that is conserved 
by H'^ + H'^). Following the Fermi's golden rule we can calculate the absorption probabihty 
for the N^{N + 1) transition: 

P{Xn^Xn+i)e = '^\H'{XnXn+i)\^S{E{Xn+i) - E{Xn) - E) (10) 
where the state |Ajv) corresponds to N Coulomb-correlated electron-hole pairs. 

IV. COUPLED CAN DOTS 

We study the multi-exciton optical response of GaN/AlN QDs in the range of parameters 
experimentally realized in [16]: the dot height is varied between 2 to 4nm and the base 
diameter is varied from 10 to 17 nm, proportionally to the height (according to experimental 
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data in [16]). The material in the dot is assumed to have a constant composition of GaN, 
while in the barrier is pure AIN, thus neglecting intermixing. The valence band discontinuity 
is set, according to experimental values, to 0.5 eV, the conduction band to 2.0 eV [17]. The 
typical system considered is composed by two QDs stacked along the growth 2;-axis (see 
Fig. 1); in the in-plane directions the confined potential is assumed to be parabolic. In 
Table I we compare the difi^erent parameters used in this paper for GaN and GaAs quantum 
dots. GaN has higher electron/hole masses and conduction/valence band discontinuities. 
The main feature of wurzite compared to zincoblende GaN heterostructures is the strong 
built-in electric field. The strength of the intrinsic field is of the same order of magnitude 
inside and outside the dot (MV/cm), but it is opposite in sign, antiparallel to the growth 
direction inside the dot. The built-in electric field in GaN QDs and AIN barriers is calculated 
according to [21]: 

T ( P'"' — P'' ^ 

_ ^bi\-^ lot ^ Lol) {'^'^) 



where ebr,(d) is the relative dielectric constant of the barrier (of the quantum dot), P^f^^^ 
is the total polarization of the barrier (of the quantum dot), and Liyr,{d) is the width of the 
barrier (the height of the dot). The value of the field in the barrier F^^ is obtained by 
exchanging the indices hr and d. Equation (11) is derived for an alternating sequence of 
quantum wells and barriers, but it is a good approximation also in the case of an array 
of similar QDs in the growth (z) direction. In our approach the in-plane components of 
the built-in electric field arc in fact "absorbed" in the strongly confining bidimensional 
parabolic potential which in addition preserves the spherical symmetry of the ground state 
[22]. Our modeling is supported by the agreement with the experimental findings in [22]. 
The polarization Pto/'^'* is the sum of the spontaneous polarization charge that accumulates 
at GaN/AlN interfaces and the piezoelectric one, Ptot — Ppiezo -^sp^*^" The piezoelectric 
charge is induced by the lattice mismatch and by the thermal strain {Pple^o — Pml*^ Pu^'^)- 
All the material parameters arc the one used in rcf [21], except for P^^ and which depend 
on the Aluminum concentration in the barrier: in the quantum wells considered in [21] such 
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concentration is Al=0.15 while for GaN/AlN quantum dots is Al=l. In particular (7|| = 
2.48%, therefore in GaN-based dots spontaneous polarization and piezoelectric fields have 
similar values while in AlQ i^GaQ^^^N / GaN quantum well the piezoelectric field contributes 
only 10% of the total value. The other parameters necessary to calculate the electric field 
are: P'^l" = -2(e33Ci3/C33 - e:,^)a^'\ Pt = -3.2 * lQ-^C/m\ = QC/m\ P^; = 
-0.081C/m2 and P^, = -0.029C/m2. 

The giant internal field strongly modifies the conduction and valence bands along the 
growth direction and causes the separation of electrons and holes, driving the first one 
toward the QD top and the latter toward its bottom. This corresponds to the creation 
of intrinsic dipoles. In this context, the charge distribution of two vertically coupled GaN 
dots occupied by one exciton each can be described as two dipoles aligned along the growth 
direction. This is exemplified in Fig.l, where we plot the electron and hole single-particle 
distributions corresponding to the lowest biexcitonic state (with parallel-spin excitons) for 
one of the GaN/AlN quantum dot nanostructure considered. Nearby quantum dots are then 
coupled by the corresponding exciton-exciton (dipole-dipole) interaction. 

V. BIEXCITONIC SHIFT AND OSCILLATOR STRENGTH 

The energy renormalization of the excitonic transition in the presence of another exciton 
is known as biexcitonic shift. A biexcitonic shift of the order of few meV in two coupled 
GaN dots is the prerequisite for the implementation of conditional quantum dynamics in the 
quantum information/computation scheme proposed in [12]. In the aforementioned scheme 
excitons in different QDs are manipulated by energy-selective addressing; additionally, due 
to the strain field, two stacked QDs are in general not identical, therefore we set the difference 
between the well widths of two stacked QDs to be 8% [24] . This allows for energy-selective 
generation of ground-state excitons in neighboring QDs. For the range of parameters con- 
sidered, the barrier width is such to prevent single-particle tunneling and to allow at the 
same time significant dipole-dipole Coulomb coupling. In contrast to GaAs quantum dots. 
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even with a 2 nm barrier the single-particle tunneling is negligible on the nanosecond time 
scale (i.e. on a time-scale comparable to the excitonic lifetime). This is due to the higher 
valence/conduction band gaps and electron/hole effective masses. 

The theoretical approach employed to study the optical response of our GaN nanostruc- 
ture is a fully three-dimensional exact-diagonalization scheme, as described in Sec. III. The 
confinement potential of GaN " macromolecule" is modeled as parabolic in the x — y plane 
and as a sequence of triangular-like potential wells along the growth [z] direction (see Fig.l). 
We evaluate the excitonic (0^1) and biexcitonic (1^2) optical spectra. The biexcitonic 
(1 — > 2) optical spectrum describes the creation of a second electron- hole pair in the presence 
of a previously generated exciton. Here, we shall consider parallel-spin configurations only. 
For all the structures considered, the two lowest optical transitions are, respectively, direct 
ground-state excitons in dot a and h (see Fig. 1). In Fig. 2 we show an example of such an 
absorption spectrum for the parameters of Fig.l. In such example a negative exciton-exciton 
couphng As = 5.7 meV is the signature of the creation of vertically aligned dipoles forming 
the biexcitonic ground state. 

Let us focus on the biexcitonic shift defined in the present case as the energy difference 
between the ground-state biexcitonic transition (given a ground-state exciton in dot a) and 
the ground-state excitonic transition of dot h. This quantity provides the essential coupling 
for realizing conditional gates in exciton-based all-optical quantum information schemes 
[10,12]. 

The biexcitonic shift can be engineered by varying the coupled GaN dots parameters 
with self- assembled growth. We analyze how it depends on the height and base diameter 
of the dot and barrier between the dots. We also study the corresponding variation of 
the oscillator strength of the transition. Figure 3 shows how the biexcitonic shift depends 
on the dimensions (height and base diameter D) of the dot, for a barrier width that 
is equal to 2.5 nm. We have considered three cross-sections of the space parameters, the 
first (solid line, A) keeping the QD height fixed {Ld = 2.5 nm), the second (dashed line, 
B) fixing the base diameter to = 10 nm (small dots), and the third (dashed-dot line, 
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C) varying the base diameter proportionally to the height, according to the experimental 
relation D = S.SL^ + 3nm [16]. The height and base diameter values plotted corresponds 
to the smaller dot. 

The excitonic dipole length is roughly proportional to the height of the dot because of the 
MV/cm built-in electric field; therefore the dependence of the exciton-exciton interaction on 
the QD height is the strongest (curve B). 

By looking at curve A, instead, we notice that the spreading of the wavefunction related 
to a wider dot basis, negatively affects the biexcitonic shift: in fact this decreases from 5.1 
to 4.3 meV, as the base diameter goes from 10 to 17 nm [16]. Exciton-exciton interaction 
is favored by "localized" states, virtually achieving a maximum in an idealized "point-like" 
particle case. 

Curve C presents the experimentally most common case in which base and height of the 
dots change simultaneously. For reahstic parameters the biexcitonic shift can be up to 20% 
smaller than in curve B, where the wavefunction is more localized, being the diameter fixed 
to the value of the "small" dots. In such case it is possible to achieve biexcitonic shifts up 
to 9.1 meV. 

Our results show that the best strategy to achieve large biexcitonic shift is to grow "high" 
and "small diameter" dots. Unfortunately the oscillator strength (OS) of the ground-state 
transition strongly decreases with the height of the dot, since it is proportional to the overlap 
of electron and hole wave functions. Since the couphng of the laser field to the considered 
transition is directly proportional to this quantity, in view of an all-optical manipulation of 
correlated electron-hole pair excitations, it is of utmost importance to study the dependence 
of the OS on the various parameters. In Fig. 4 the OS corresponding to the excitonic ground 
state of dot h is plotted for the same parameters of Fig. 3 (same labeling of the curves). 
Curve B (fixed base) shows that the OS decreases super-exponentially with the height of the 
dot. It changes over three order of magnitude from 2 nm to 4 nm height dot. It is interesting 
to notice that, on the other hand, the width of the dot does not influence the electron-hole 
overlap, so curve A (fixed Ld) is practically a constant over the whole D range and the two 

13 



curve (B) and (C) corresponding to different diameter, coincide [25] . 

In the range of height values considered in Fig. 3, the OS varies over three orders of 
magnitude, so care must be taken in a future quantum information processing experiment 
in order to optimize at the same time biexcitonic shift and OS [26] . 

Let us now examine the influence of the barrier on the biexcitonic shift and oscillator 
strength. In fig. 5 we show (upper panel, curve A) the decreasing of the biexcitonic shift 
with the barrier width for two coupled dots of 2.5 nm and 2.7 nm of height. We notice that it 
deviates from a dipole-dipole interaction behavior (curve B) and it is 0,5 meV even at 10 nm 
barrier. The reason of this behavior can be inferred from equation (1). The electric field in 
the dot increases with the barrier width, and, as a consequence, also the dipole length. The 
two competing effects result in a slower decreasing of the biexcitonic shift with the barrier. 
Curve B corresponds to the point- like approximation, for a fixed dipole length (about 1.5 
nm) corresponding to a 2 nm barrier. We see that the "point- like" biexcitonic shift is two 
times bigger when the barrier is 2 nm, there is a crossover at 4 nm, after which it remains 
smaller than the exact result. The most striking consequence of the barrier- dependent huilt- 
in electric field, is that the oscillator strength strongly varies too (lower panel), even if the 
dimensions of the two dots are kept fixed. In addition (and in contrast to Fig. 3 and 4), in 
this case both OS and biexcitonic shift are increased by a smaller barrier, suggesting that the 
best structure for optical excitonic manipulations and energy-selective addressing techniques 
is characterized by small inter-dot barriers. In this respect it would be very interesting to 
have experimental confirmation of our prediction. 

In Fig. 6 we plot the direct comparison of electron and hole distributions of the 2.7 nm 
dot, for two different barrier width (2 nm and 10 nm). The width of the barrier changes 
the internal electric field in the dot, that in turns strongly modifies the wavefunctions. The 
dipole length changes of about 30% (Fig. 6), going from about 1.5 nm at 2 nm barrier, to 
1.9 nm when the barrier is 10 nm wide. 

Compared to GaAs quantum dots, nitride dots have higher effective masses m* of both 
electrons and holes, lower dielectric constant and higher (about two times) conduction and 
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valence band discontinuities Ve^h- As a consequence, the GaN bulk excitonic Bohr radius 
is about four times smaller than the GaAs one and the distributions of confined particles 
are more "localized", because their wavelength penetration L^h, where L^h = i ^ =, 
depends on such parameters. 

The reduction of these two characteristic lengths causes a faster decrease of the electron- 
hole overlap, as a function of the dipole length and consequently a stronger sensitivity 
of the oscillator strength (roughly proportional to the electron-hole overlap in the strong 
confinement regime we are considering) to quantities such as the height of the dot (Fig. 4) 
or the different built-in electric field (Fig. 5), which strongly affect the dipole length. 

Our results demonstrate that there exist a wide range of parameters for which the biex- 
citonic shift of the order of meV. This is a central prerequisite for realizing energy-selective 
addressing with sub-picosecond laser pulses, as requested, for example, by all-optical quan- 
tum information processing schemes [10,12] or read-out devices [13]. 

VI. BLUESHIFT AND ABSENCE OF MOLECULAR STATES 

Another interesting effect peculiar of hexagonal GaN/AlN quantum dot is the blueshift 
of the ground state transition when the distance between the dots is decreased, without the 
hfting of the degeneracy of bonding-antibonding states. In GaAs-based quantum dots [29] 
there is a redshift and an increasing energy difference between bonding and anti-bonding 
states, which are spread over the whole macromolecule for both electron and hole; on the 
contrary in GaN QDs, over the range of parameters here considered, the lowest states pre- 
serve their atomic- like shape. This depends on the fact that both electron and hole effective 
masses and valence/conduction-band discontinuities are much larger than in GaAs, there- 
fore decreasing the atomic-like wavef unction overlap responsible for the molecular bonding. 
However, as discussed in the previous section, the built-in electric field is not only a function 
of the dot parameters, but it depends strongly on the barrier width. In particular it increases 
when the barrier is increased, to saturate at the value ~ {P^^ — Pf^^) / (eoe^) , corresponding 
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to an isolated dot. For examples it increases from 5 MV/cm for a 2 nm barrier to 9 MV/cm 
for a 10 nm barrier for the dot parameters of Fig. 6 (see inset). In Fig. 7 we consider two 
identical dots separated by increasing barrier width. As a consequence of the electric field 
change, the ground state energy changes of more than 600 meV, from 3337,5 meV to 2709,5 
meV (upper panel of Fig. 7). Also the oscillator strength is modified (lower panel of fig. 7) 
of one order of magnitude, decreasing very fast up to 5 nm barrier width, showing a slower 
decreasing afterward. In the particular case of exactly identical quantum dots, because of 
the symmetry of the system, the single exciton eigenstates are {ipa) = "^(|01) + |10)) and 
= :^(|01) — |10)) and the energy splitting is given by 2Vf where Vp is the Foerster 
energy [30] . Vp is proportional to the square of electron-hole overlap [30] , therefore is much 
lower in GaN than in GaAs. For example the maximum value for the parameter considered 
in this work, corresponding to 2 nm dot's height and 2 nm barrier, is only 0.05-0.06 meV. 
However we stress that when the dots are slightly different, this effect rapidly vanishes since 
it depends on the ratio (Vp/A)^ where A is the energy difference between the two levels 
coupled by the Foerster transfer process. In GaN dots, differently form GaAs, with an 8% 
difference in size, A is already of the order of hundreds of meV, causing this effect to be 
extremely difficult (if not impossible) to be observed experimentally. 

To qualitatively investigate the effects of the different parameters involved we have 
calculated the ground state electron and hole distribution for some GaAs-GaN "mixed 
case", i.e. some artificially designed molecules with one GaN parameter (electric field, 
valence/conduction band offsets, effective masses) substituted by the corresponding GaAs 
one. We will consider a nanostructure composed by two dots of 2.5 nm and 2.7 nm separated 
by a 2 nm barrier. As a benchmark we also plot the "pure GaAs" case (fig. 8b): its ground 
state has both hole and electron distributions delocalized over the macromolecule, forming 
bonding molecular states. The corresponding "pure GaN" case is plotted in Fig.l and it 
is characterized by electron and hole distributions localized in a single dot. The coupled 
GaN dots without the giant electric field (fig.8d ) still present an atomic character. The 
same holds for an artificial GaN macromolecule with GaAs effective masses (Fig.8c ), and 
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an artificial GaN nanostructure with GaAs conduction/valence band structure (Fig.Sa). It 
is interesting to note that the electron and hole distribution of a GaN dots without electric 
field (Fig. 8d) are more extended than when a giant electric field is present (Fig.l). We con- 
clude that the single dot confinement of the excitonic wavefunction in GaN QDs is a robust 
feature and that all the parameters involved are responsible for the absence of molecular-like 
states in GaN coupled dots. 

VII. SUMMARY AND CONCLUSIONS 

We have performed a detailed investigation of the optical spectrum of coupled GaN 
quantum dots. In particular we have shown some effects peculiar to wurzite nitride materials 
compared to GaAs-based nanostructures: the absence of ground state exciton wavefunction 
delocahzation even for relatively short barriers (up to 2 nm) due to the large effective 
mass/band offsets of nitrides; the presence of a large exciton-exciton interaction between 
neighbor quantum dots, caused by the giant intrinsic electric fields. We have also shown 
that a shift of energy levels of identical quantum dots is expected when the distance between 
them is varied: such blueshift is not caused by molecular coupling of the wavefunction, but 
by the decreasing of the built in electric field as the inter-dot barrier decreases. We have also 
shown how it is possible to engineer the inter-dot biexcitonic shift and the corresponding 
oscillator strength by varying the structural parameters (base, height, barrier) of the dots. 
Such analysis is crucial in the mainframe of quantum information processing schemes [12] 
and in general for all-optical devices based on energy-selective addressing [13] in which the 
conditional excitonic dynamics is based on such quantity. 
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FIGURES 

FIG. 1. Electron (upper panel) and hole (lower panel) particle distribution (dotted line), con- 

dTiction (upper panel) and valence (lower panel) band structure (solid line) along the growth 
direction for two coupled GaN dots of, respectively, 2.5 nm and 2.7 nm of height, separated by a 
2 nm AIN barrier. 

FIG. 2. Excitonic (solid line) and biexcitonic (dashed line) absorption spectrum of the GaN 
coupled dots of fig.l. 

FIG. 3. Biexcitonic shift of the ground state transition in dot b for two coupled GaN dots as a 
function of dot height and base diameter. In curve (A) only the base diameter of the dots is changed 
{D = lOnm); in curve (B) only the height of the dots is changed {L^ = 2.5nm), while in curve 
(C) D is varied proportionally to L^. The parameters for the parabolic potential is varied in the 
range: htOg = 74 290meV, htOh = 33 130meV. In particular the values analized in the figures 
are D=10nm {TiLOe = 290meV,hLOh = ISOmeV), D=11.75nm(^a;e = llOmeVfiLOh = 76meV), 
D=13.5nm(^a;e = ISlmeV ,TiLOh = 59meV), D= 15.25nm(^a;e = 97meV,TiLOh = ASmeV), 
D=17nm(?iu;e = T^meV fiujh = 2>2>meV) 

FIG. 4. Oscillator strength of the ground state transition in dot b for the same parameters of 
Fig. 3. Labelling as in Fig. 3. 

FIG. 5. Upper panel: curve (A) shows the biexcitonic shift of the ground state transition in dot 
b for two coupled GaN dots of 2.5 and 2.7 nm of height vs barrier width between the dots. Curve 
(B) shows the corresponding biexcitonic shift in the point-like charge approximation. Inset: as in 
main panel, but for large barrier widths. Lower panel: oscillator strength for the same parameters 
described above. Inset: internal electric field in dot b vs barrier width. 

FIG. 6. Excitonic distribution in a 2.5 nm dot of a coupled dot nanostructure when the barrier 
width is 2 nm (solid line) and 10 nm (dotted line). 
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FIG. 7. Exciton ground state energy (upper panel) and oscillator strength (lower panel) in the 
dot b of two identical coupled GaN dots of 2.5 nm of height vs barrier between the dots. 

FIG. 8. Electron and hole particle distribution along the growth direction for two coupled dots 

of, respectively, 2.5 nm and 2.7 nm of height, separated by a 2.5 nm AIN barrier. The dot parameter 
sets corresponds to: a)GaN, except valence/conduction band offset (GaAs one); b)GaAs; c)GaN, 
except electron/hole masses (GaAs ones); d)GaN, except that the electric field is set to zero. 
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TABLES 



GaN GaAs 

EgapieV) 3.5 1.519 

m*(mo units) 0.2 0.067 

ml^{mo units) 1 0.34 

AEc{eY) 2 1 

AEy{eV) 0.5 0.3 

e(relative) 9.6 12.10 



TABLE I. Material parameters used in the calculations 
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Fig.l S. De Rinaldis et al. 
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